Field-induced negative differential spin lifetime in silicon 
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We show that the electric field-induced thermal asymmetry between the electron and lattice 
systems in pure silicon substantially impacts the identity of the dominant spin relaxation mechanism. 
Comparison of empirical results from long-distance spin transport devices with detailed Monte- 
Carlo simulations confirms a strong spin depolarization beyond what is expected from the standard 
Elliott- Yafet theory already at low temperatures. The enhanced spin-flip mechanism is attributed 
to phonon emission processes during which electrons are scattered between conduction band valleys 
that reside on different crystal axes. This leads to anomalous behavior, where (beyond a critical 
field) reduction of the transit time between spin-injector and spin-detector is accompanied by a 
counterintuitive reduction in spin polarization and an apparent negative spin lifetime. 



In compound semiconductors, the eventual reduction 
in drift velocity of conduction electrons with increasing 
applied electric field is known as negative differential mo- 
bility or the Gunn effect [l|, 0] . In this field regime (typi- 
cally several kV/cm), hot electrons scatter into low- lying 
secondary energy minima in the conduction band where 
the effective mass is larger, reducing their kinetic energy. 
The multivalley band structure of silicon also allows for 
the existence of this phenomenon but only at low temper- 
atures; for all T >30 K, the drift velocity increases and 
eventually saturates with increasing applied field [HjH- 
Therefore, at elevated temperatures the time-of-flight of 
conduction electrons across the Si channel of a transport 
device drops monotonically with increasing electric field. 
If electrons are initially spin polarized, then the accepted 
Elliott- Yafet spin relaxation theory suggests that the spin 
depolarization during transport is dependent only on the 
time-of-flight. In this theory the spin and momentum 
relaxation times are proportional [a, @, so the result- 
ing spin polarization increases with electron drift veloc- 
ity. Indeed, we have confirmed this expectation in previ- 
ous experiments where low and moderate applied fields 
(<1 kV/cm) were used in studying the extraordinarily 
long spin lifetime [3,l8|, and in demonstrating spin injec- 
tion and detection in ferromagnet-silicon hybrid systems 

In this Letter, we experimentally demonstrate an unex- 
pected dependence of the spin polarization on the electric 
field in silicon at T >30 K in high electric fields. With 
increasing field, the spin polarization of detected elec- 
trons first increases as expected from the Elliott- Yafet 
static lifetime model; however, above 2 kV/cm it starts 
to decrease, showing a Gunn-effect dependence akin to 
a negative differential spin lifetime without any simul- 
taneous negative differential charge/spin mobility. The 
origin of this counterintuitive behavior is then elucidated 
by Monte Carlo simulations and a quantitative analytical 
description. When the electron ensemble is out of ther- 
mal equilibrium with the lattice, an efficient spin relax- 



ation mechanism becomes accessible due to field-induced 
intervalley scattering. We quantify the spin relaxation 
time as a function of both the lattice and electron en- 
semble temperature. The latter provides a means to de- 
termine the dependence of spin relaxation in silicon on 
the electric field and will enable the optimization of spin- 
tronics devices. 

Coherent spin precession and spin valve measurements 
were performed to observe the nonequilibrium depolar- 
ization effect and to quantify the negative differential 
region of spin lifetime. In both experiments, we em- 
ployed all-electrical devices in which spin-polarized elec- 
trons (aligned with the in-plane magnetization direction 
of a ferromagnetic thin-film source) are tunnel injected 
through a Schottky metal contact and into a 225 /im - 
thick wafer of nominally undoped Si(100). The electrons 
then drift across the wafer thickness due to an applied 
electric field, and are collected by a second ferromag- 
netic film where their spin is analyzed using a ballistic 
spin detection scheme. The results presented below do 
not depend on the injection and detection techniques but 
only on the spin and charge transport characteristics of 
the Si channel. We therefore include all device-specific 
description in the supplemental material Uj, and refer 
the interested reader to Refs. 0,SII3| for further details. 

The time-of-flight distribution of the electron current 
can be recovered from quasistatic spin precession mea- 
surements by applying an external magnetic field, B, 
perpendicular to the injected spin direction but paral- 
lel to the electric field [121. This magnetic field induces 
spin precession at frequency lo — gfiBB/h, where g is the 
electron g- factor, {1b is the Bohr magneton, and Ti is the 
reduced Planck constant. We denote the time-of-flight 
distribution by D(t) where its mean and standard devi- 
ation are, respectively, measures of the average transit 
time and of diffusion and dephasing effects in the chan- 
nel. The signal contribution from the spin component 
parallel to the detector magnetization of an electron ar- 
riving at the detector in the time interval [i, t + dt] is 
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FIG. 1. Average transit time across a 225-micron-thick in- 
trinsic silicon device as a function of the applied voltage for 
various temperatures. Error bars indicate the transit time 
uncertainty (extracted from the width of the time-of-fiight 
distribution; see lower inset). Top inset: Symmetrized spin 
precession data at T = 61 K and 80 V (3.5 kV/cm). The 
data show high spin coherence with well-defined oscillation 
field period. Lower inset: spin current transit time distribu- 
tion obtained by transforming the precession signal (see text). 
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FIG. 2. Measured current polarization in a ferromag- 
net/silicon/ferromagnet device as a function of transit time 
in the silicon channel. Solid lines indicate exponential fits 
to the low-field (long transit time) data and indicate long 
spin lifetimes in that regime. Effects of spin depolarization 
from electric-field-induced spin relaxation are evident at high- 
fields (short transit time; circled data). Inset: Example spec- 
troscopy at 61 K showing spin signal in a parallel (P) and 
antiparallel (AP) magnetic configuration, with background 
signal for subtraction. 



therefore D(t) cos uitdt. The variation in quasistatic de- 
tected signal is then D(u>) oc J Q D(t) cos uttdt; by repeat- 
ing the measurement at various applied magnetic fields 
one can map theprecession frequency dependence of the 
detected signal [8j. Finally, the empirical time-of- flight 
distribution is recovered without any model fitting by the 
inverse Fourier transform of D(u>). An example of this 
transformation between D(uj) and D(t) is shown in the 
coupled insets to Fig. [TJ The main figure shows the aver- 
age transit time across the silicon channel as a function 
of applied voltage for several temperatures. Clearly, in- 
creasing the internal electric field with applied voltage 
reduces the transit time until the onset of velocity satu- 
ration for voltages > 60 V (electric field ~ 2.7 kV/cm) 
It is important to note that no negative differential 
mobility is seen. 

To measure the spin polarization of collected electrons, 
we have performed independent spin- valve measurements 
as a function of the electric field. Here, a small external 
magnetic field of w 20 Oe is applied along the source mag- 
netization axis and thus no spin precession is induced. 
The final spin polarization after transport is extracted 
by the ratio V = (Ip — Iap)/(Ip + Iap), where Ip is 
the measured signal current in a configuration where the 
in-plane injector and detector magnetization directions 
(and hence spin initialization and measurement axes) are 
parallel, and Iap is for antiparallel configuration. The 
inset in Fig. 2(a) shows an example of this spectroscopy 
taken by interleaving P and AP measurements at each 
applied voltage to avoid signal drift from field-induced 



stress. We also include the background detector current 
taken under conditions of zero injection current after the 
signal measurement. It has subsequently been subtracted 
in the polarization calculation to avoid misinterpreting a 
spurious dilution for spin depolarization. 

Figure HJa) shows the measured polarization as a func- 
tion of the average transit time (T tr ). This ratio de- 
pends on the spin relaxation time in the Si channel by 
(I p -Iap)/(Ip + Iap) = V a er Tt -l T " whereto (limited by 
the spin-injection and detection efficiencies of the device) 
is the optimal attainable value. The figure shows that at 
long transit times, the polarization increases with reduc- 
ing the transit time, as expected from the Elliott- Yafet 
theory. However, at short transit times (circled data) 
the trend is unexpectedly opposite. This observation of 
a nonmonotonic spin polarization Gunn effect is the main 
experimental result of this Letter. 

The origin of this phenomenon is a transition to a 
previously-ignored regime where electric field directly en- 
ables a spin relaxation pathway. The field-induced mo- 
mentum relaxation enhancement, as implied by the satu- 
ration in charge transport data of Fig. 1, is not commen- 
surate with the spin relaxation enhancement. Applying 
the accepted Elliott- Yafet theory (proportionality of spin 
and momentum relaxation times) would therefore lead to 
the false conclusion that the rising polarization with ini- 
tially increasing transit time is indicative of an unphysical 
negative spin lifetime. 

We have performed Monte-Carlo simulations in order 
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FIG. 3. (a) Electron drift velocity, (b) mean energy, (c) spin 
relaxation time and (d) final polarization as a function of elec- 
tric field, calculated from numerical integration of the distri- 
bution obtained from Monte Carlo simulation. The dotted 
lines in (c) denote Eq. (1). 



to elucidate the charge transport and spin relaxation of 
conduction electrons heated by the electric field ("hot" 
electrons). A full description of the numerical procedure 
is provided in the supplemental material 11], and here 
we summarize the important features. Ellipsoidal energy 
bands are used to model the equivalent six conduction 
band valleys fl3| . Momentum relaxation mechanisms 
are modeled by electron-phonon interactions (both in- 
travalley and intervalley processes) and intravalley elec- 
tron scattering from ionized impurities 14]. Between 
scattering events, electrons are treated as classical parti- 
cles accelerated by the electric field. Typically, an out- 
of-equilibrium electron distribution reaches its steady- 
state within 1 ns regardless of the initial condition. Fig- 
ure EJa)-(b) show the corresponding drift velocity and 
mean energy as a function of applied electric field. Hot 
(cold) valleys refer to the four (two) valleys whose axis is 
perpendicular to (collinear with) the electric field. The 
mean energy in a hot (cold) valley is higher (lower) due to 
the different projections of electric field on the ellipsoidal 
energy bands. 

The spin relaxation due to electron-phonon interac- 
tions is calculated by integration of intravalley and in- 
tervalley spin-flip matrix elements [lij while using the 
Monte-Carlo hot-electron distributions. The solid lines 
in Fig. |3][c) show the results of this numerical proce- 
dure. Note that spin relaxation due to scattering with 
ionized impurities is negligible in nearly intrinsic wafers. 
As a result the total spin lifetime in our devices is signif- 
icantly longer than in heavily doped Si channels [1 6l — ll 8j j . 
FigureEJd) shows the spin polarization Vq exp(— r tr /r s ), 
where Vq = 0.125 is chosen to fit the experimental injec- 
tion and detection efficiencies. The average transit time 
r tr is calculated from the drift velocity after transport 



across 225 ijiii. At low fields, the polarization rises with 
electric field since the increase of drift velocity surpasses 
the decrease of spin relaxation time. As the drift ve- 
locity begins to saturate in high fields, the polarization 
drops slowly due to the enhanced reduction of the spin 
relaxation time. This dependence of spin polarization on 
the electric field agrees well with the experimental results 
and reproduces the Gunn-type behavior (here shown as 
a function of the field). 

We focus on the important underlying spin relaxation 
mechanism and analytically quantify the observed ef- 
fect. From the mean energy (Fig. |3jb)) , one can see that 
electrons driven by the electric field become hot enough 
to undergo intervalley electron-phonon processes during 
which the electron delivers to the lattice a few tens of 
meV We consider /-processes at which electrons are 
scattered between valleys of different crystal axes. This 
process dominates the spin relaxation of hot electrons 
since it involves a direct coupling of valence and con- 
duction bands [15|. To conserve crystal momentum, the 
phonon wavevector resides on the S axis. The symmetry- 
allowed phonon modes for spin relaxation are Si and S3 
with respective phonon energies of « 47 meV and 
f2/,3 ~ 23 meV. The S3 mode allows for scattering be- 
tween all valleys and the Si mode restricts them to the 
case that one of the involved valley axes is not perpendic- 
ular to the spin quantization axis. To analytically quan- 
tify the spin relaxation we functionalize the hot electron 
distributions. Fig. HJa) shows the Monte-Carlo steady- 
state energy distributions in hot and cold valleys at 30 K 
and 4 kV/cm. The distribution (in each of the valleys) 
can be described by a two-component heated Boltzmann 
distribution. At the low energy part, the effective tem- 
perature of the electron distribution can be extracted 
from the mean energy |/cB?e, shown in Fig. 0(b) . At 
the higher energy part, intervalley process tend to cool 
the system. To simplify the analysis below we employ 
an effective electron temperature, T' e = T + 7(T e — T) 
where T is the lattice temperature and 7 ~ 0.9 is a con- 
stant that mimics the cooling effect due to intervalley 
scattering at high electron energies [ll[ . We denote effec- 
tive lattice and electron parameters by yi = fi/^/fcsT, 
y'i „ = flf t i/kBT{, ^ where i denotes the phonon modes 
and T' e is the effective temperature of the electrons in 
a cold (/1 = c) or a hot (/1 = h) valley. Using the above, 
we arrive at an analytical spin lifetime (ll| . 
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where C = 0.036 ns 1 is a constant related to the spin- 
orbit coupling parameter of the X point at the edge 
of the Brillouin zone. Ai, h = 8(12), A 3 , h = 1.5(1.25), 
A\,c = 8 (4) and Az. c = 0.5 (0.75) are symmetry related 
parameters when the electric field is collinear with (per- 
pendicular to) the spin-quantization axis. The n c and 
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FIG. 4. (a) Electron distributions in hot (solid line) and cold 
(dashed line) valleys. The total electron density is 10 12 cm -2 , 
the electric field is 4 kV/cm and the lattice temperature is 
30 K. (b) Ratio between electron densities in cold and hot 
valleys as a function of the field, (c) Experimental depolar- 
ization at high-fields (extracted from Fig. (2)). (d) Charac- 
teristic scale of the drop in spin polarization as a function of 
temperature. 



rih denote, respectively, the fractional population at cold 
and hot valleys where 2n c + Arih = 1. Figure SJb) shows 
the repopulation ratio, n c /n/j. The asymmetry in val- 
ley population is largest at intermediate fields since elec- 
trons that reside in hot valleys become energetic enough 
for intervalley scattering to cold valleys. At high fields, 
scattering in the opposite direction also becomes ac- 
cessible and the valley population is more symmetrical 
{fic/rih 2 in this regime). Substituting the extracted 
values of and T' in Eq. (JT]) reproduces the spin re- 
laxation of hot electrons as can be seen from the compar- 
ison between the dotted and solid lines in Fig. ^c). At 
equilibrium conditions where T' e ^ ~T this mechanism is 
greatly suppressed (especially at low temperatures). 

We can quantitatively compare the results of our cal- 
culations with empirical data by extracting a character- 
istic voltage scale Vo for /-process-induced spin depolar- 
ization, where we approximate V w Vo(l — ^r) in the 
high electric-field regime. Fits to the measured spin po- 
larization data are shown in Fig. |4tc), and the consis- 
tent temperature dependence of Vo for several devices 
is shown in Fig. 0Jd). The Monte-Carlo prediction ex- 
tracted from the high-field regime in Fig. ISJd) closely 
resembles the empirical values in both magnitude and 
lattice-temperature dependence due to the more efficient 
generation of intervalley scattering (and hence lower Vo) 
at lower temperatures. By taking the high-field limit of 
Eq. (dJ, one can write V « v d /[d(T- 1 )/dE] w 0.7 kV 
where Vd is the saturated drift velocity. This close corre- 
spondence confirms our interpretation of field-induced /- 



process spin depolarization in the experiment, in a regime 
where acoustic-phonon-mediated scattering as well as 
scattering with states at the spin hot-spot [HI, [l9| are 
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too small to account for this effect 

Finally, we note that a recent theoretical proposal sug- 
gesting that stochastic polarization fluctuations can be 
amplified by spin-dependent mobility [ioj] has also been 
termed a "spin Gunn effect" . Our experiment and the- 
ory differ from this scheme in that the mobility and dif- 
fusion constants are spin independent, electron-electron 
collisions are negligible [21], and the physical origin of 
the effect is attributed to the signature of spin-orbit cou- 
pling on electron-phonon intervalley scattering. The phe- 
nomenon observed here, as in its charge-based counter- 
part, is due to a strong electric field- induced relaxation 
which leads to a qualitatively different spin transport 
regime distinct from expectations based on the Elliott- 
Yafet theory. 

In conclusion, high electric fields present in silicon 
devices can substantially change the dominant physi- 
cal mechanism of spin relaxation. In this regime, the 
Elliott- Yafet mechanism mediated by intravalley acous- 
tic phonons is far outweighed by the depolarizing effects 
of inelastic scattering with intervalley /-process phonons 
created by the efforts of the system to recover thermal 
equilibrium. This behavior is expected to be critically 
important in the design of devices making use of spins 
to transmit information, especially when strong static 
electric fields are required [22j . Similarly, the derived de- 
pendence of the spin lifetime on the electric field is of 
fundamental importance to the design of semiconductor 
devices that make use of spin as an alternative degree of 
freedom [23l426l | . 
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